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We report results of the first phase of observations with Australia Telescope Com- 
pact Array (ATCA) at 5 and 9 GHz of the fields around 411 7-ray sources with 
declinations < +10° detected with Fermi but marked as unassociated in the 2FGL 
catalogue. Wc have detected 416 sources with flux densities in a range of 2 mJy to 
6 Jy in the fields of 267 7-ray sources within their position error ellipses drawn to 
cover the area of 95 per cent probability of their localisation. Of them, 142 objects 
were detected in both bands. We found 95 flat spectrum sources in our sample. The 
majority of them are weak and for this reason were not found in previous surveys. 
Approximately 1/4 of our sample, 107 objects, have the probability of being associ- 
ated by more than 10 times than the probability of being a background source found 
in the vicinity of a 7-ray object by chance. We present the catalogue of positions of 
these sources, estimates of their flux densities and spectral indices when available. 
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1 INTRODUCTION 

Analysis of two years of Fermi/LAT observations yielded a 
catalogue of 1872 7-ray sources (Abdo et al. 2012, thereafter 
2FGL). Of them, 70 per cent have associations with blazars, 
pulsars, supernova remnants, and other objects. The proce- 
dure for assigning associations is described in full detail in 
Ackermann et al. (2011). However, 2FGL catalogue still did 
not provide associations for 573 objects. Position accuracy 
of Fermi ranges from 0.1 arcmin to 16.6 arcmin with the me- 
dian la uncertainty 2.0 arcmin. These large position errors 
prevent finding associations by direct matching the 2FGL 
against optical or infrared catalogues. 

Previous analysis of IFGL catalogue (Kovalev 2009) 
has confirmed earlier EGRET results: 7-ray emission and 
parsec-scale radio emission are strongly related. Extending 
this study to the 2FGL catalogue, we found that 745 out of 
1872 Fermi sources, roughly one half, have been detected in 
VLBI surveys at 8 GHz by December 2012. All these objects 
are active galactic nuclei (AGN). This is the major popula- 
tion of point-like Fermi sources outside of the Galactic plane. 
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Other type of objects associated with 7-ray sources are su- 
pernova remnants, novae, pulsar wind nebulae, x-ray bina- 
ries, microquasars, and pulsars. Unassociated sources may 
belong to any of these classes or may in part constitute an 
as-yet-unknown population. We find that high resolution ra- 
dio observation are helpful for classification of unassociated 
sources. Supernova remnants and pulsar wind nebulae are 
extended objects and high resolution radio interferometric 
observations tend to resolve out their emission making them 
undetectable. Pulsars are generally weak at high frequencies 
(> 5 GHz). Therefore, sources detected with VLBI that are 
brighter than 1 mJy are almost always AGNs. Detection of 
a radio bright AGN within an error ellipse of Fermi posi- 
tion presents a strong argument that this is likely the same 
object. Therefore, systematic high resolution radio observa- 
tions of unassociated Fermi objects promise to find all 7-ray 
sources associated with radio-loud AGNs and shrink the list 
of objects that are still not associated. 

It was shown by Lister et al. (2009) and Kovalev (2009) 
that 7-ray fluxes and radio flux densities at 8 GHz from re- 
gions smaller than 5 mas are correlated. Source variability 
conflrms that both radio and 7-ray emission comes from 
regions of parsec-scale. Therefore, a flux-limited sample of 
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Figure 1. Dependence of the decimal logarithm of the number of 
compact sources versus the decimal logarithm of their minimum 
correlated flux density from regions smaller than 50 mas at 8 GHz. 
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Figure 2. Maximum flux density at 8 GHz of a background com- 
pact source to be found with the probability of 10 per cent within 
an error ellipse with a given la standard deviation. 



7-ray AGNs, such as 2FGL, should be contained within a 
flux-limited sample of compact radio sources. Compact ra- 
dio sources from the flux-limited catalogue found within er- 
ror ellipses of sources from a 7-ray flux limited catalogue 
can be considered with a high confidence as the same ob- 
jects, provided the probability of finding a background radio- 
loud source within the error ellipse is small. The smaller the 
search area, the weaker radio-loud AGN can be associated 
with a 7-ray source. 

In order to evaluate the probability of finding a 
background source in a given search area, we investi- 
gated the cumulative all-sky catalogue of compact ra- 
dio sources detected with VLBI in the absolute as- 
trometry mode (Petrov and Kovalev, in preparation). 
This catalogue^ has 7216 objects detected in numerous 
VLBI surveys over the last several decades; Very Long 
Baseline Array (VLBA) Calibrator Survey (Beasley et al. 
2002; Fomalont et al. 2003; Petrov et al. 2005, 2006; 
Kovalev et al. 2007; Petrov et al. 2008), Long Baseline Ar- 
ray (LBA) cahbrator survey (Petrov et al. 2011b), VLBA 
Galactic Plane Survey (Petrov et al. 2011a), European 
VLBI Network (EVN) Galactic plane survey (Petrov 2012), 
VLBA imaging and polarimetry survey (Taylor et al. 2007; 
Petrov & Taylor 2011), regular VLBA geodetic observations 
(Petrov et al. 2009; Pushkarev & Kovalev 2012), and ongo- 
ing VLBI observations of 2MASS galaxies (Condon et al. 
2011). Figure 1 shows the dependence of the logarithm of a 
number of sources at all sky with flux density greater than 
5 as a function of logarithm of S determined as the median 
correlated flux density at baseline projection lengths in a 
range of 100-900 km at 8 GHz. We see that the dependence 
can be approximated as N{S) ~ 327 • S~^-'^^''' in a range of 
[0.18, 5] Jy. We interpret the deviation of log A''(log 5) from 
the straight line below 180 mjy as evidence of incomplete- 
ness. Assuming the parent population remains the same for 
sources as weak as 1 mJy, we can extrapolate the dependence 
of the number of compact sources at the celestial sphere de- 
rived for the range [0.18, 5] Jy to the range [1, 180] mJy. 

Let us consider a search for a radio counterpart within 
the area where the probability of source localisation is P. 

^ Available at http://astrogeo.org/rfc 



Then, assuming localisation errors follow the 2D normal dis- 
tribution with the second moment a, the mathematical ex- 
pectation of the number of background sources that can be 
found in that area is 

2 

M{a,S,P) = -^\n{l-P)N{S). (1) 

ZTT 

The function M{a, S) depends on both the position un- 
certainty, the flux density, and the probability of localisa- 
tion. When M -C 1, it can be interpreted as the probability 
that at least one source will be found within the area of po- 
sition uncertainty. If we fix M{a, S) to a specific value, for 
instance, 0.1, and fix localisation probability, f.e. to 0.95, we 
can find the dependence of the maximum flux density S of 
a background source that can be found within the area of 
localisation with a certain probability (10 per cent in our 
example) on the standard deviation of localisation. This de- 
pendence is shown in Figure 2. For the median la position 
uncertainty of unassociated sources from the 2FGL cata- 
logue 3.0 arcmin, the probability to find a background source 
brighter than 10.9 mJy is 10 per cent. For 80 per cents of 
2FGL source, la position uncertainty is less than 4.0 arcmin. 
The fiux density that corresponds to 10 per cent probability 
to find a background source with such a position uncertainty 
is 17.3 mJy. 

In order to check the validity of reported 2FGL po- 
sition uncertainties, we computed arcs between 745 7-ray 
objects that have counterparts with radio sources observed 
with VLBI and normalised them to their standard devia- 
tions derived from parameters of their 2FGL error ellipses: 



^- -Jal^^sin^P + al-^^cos^P 

n = ^A'2acos^5 + A^S-^ , (2) 

amaj amin 

where angle /3 is 
AS 

/3 = arctan (7r/2 - 6*) (3) 

Aa cos 5 

and CTmaj, Cmin, and 9 are semi-major, semi-minor axes and 
position angle of the 2D Normal distribution that describes 
errors of source localisation. We derived amaj, cTmin from 
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Figure 3. Distribution of position offsets normalised by their 
standard deviations between 745 Fermi sources and their coun- 
terparts detected with VLBI. The solid blue lines shows Rayleigh 
distribution with tr = 1. 



the so-called "95% confidence source semi-major and semi- 
minor axes" reported in table 3 of 2FGL catalogue by scaling 
them by 1/e. The distribution of normalised arcs between 
VLBI and 2FGL positions follows very closely Rayleigh dis- 
tribution with (7 = 1 (Figure 3). The histogram is best fit 
into Rayleigh distribution with a = 1.026. Since VLBI posi- 
tions are five orders of magnitude more accurate than posi- 
tion reported in 2FGL, they can be considered as true. We 
conclude that 2FGL position errors are realistic and can be 
used for a probabilistic inference. 

These estimates show that radio observations are an ef- 
ficient tool to identify AGNs associated with 7-ray sources. 
Nowadays, it is not difficult to detect sources with fiux den- 
sities of several mjy using radio interferometers. As we have 
shown, the cumulative catalogue of compact radio sources 
detected with VLBI is complete only to the level of 180 mJy. 
Sources with the correlated fiux densities 2-180 mJy may 
not have been marked as associations with 2FGL objects 
because they were missing from VLBI surveys. Not all as- 
sociations of 7-ray sources in 2FGL are made on a basis of 
VLBI observations. In the absence of VLBI observations, 
the source spectrum can be used as a proxy. It was found a 
long time ago that sources with spectral index fiatter than 
—0.5 tend to be compact with synchrotron radiation dom- 
inated by emission from the core. The problem is that ex- 
isting catalogues are not complete enough to determine the 
spectrum in a range of flux densities 1-100 mJy, especially 
in the southern hemisphere. 

These considerations prompted us to propose a program 
of observations of all unassociated 2FGL sources. The even- 
tual goal of the program is to find all AGNs with correlated 
fiux densities brighter than 2 mJy at 8 GHz from regions 
smaller than 50 mas within the areas of 95 per cent proba- 
bility of Fermi localisation uncertainties. Association with a 
VLBI source automatically improves positions of 7-ray ob- 
jects from arcminutes to milliarcseconds, which allows for an 
unambiguous association with sources at other wavelengths, 
for instance, optical. Combined with results of on-going pro- 
grams to find pulsars associated with Fermi sources, for in- 
stance, Barr et al. (2013). we expect the list of unassociated 
sources will shrink. At the moment we do not know whether 




Figure 4. The distribution of 573 2FGL unassociated sources. 
The area of this survey is shown with grey colour. 



the list will shrink to zero or a population of radio-quiet 
7-ray sources will be found. 

In this paper we report early results of the first step of 
the program of observations of Fermi unassociated sources: 
observations with the Australia Telescope Compact Array 
(ATCA) of 411 sources with 5 < +10°. A brief overview of 
the program is given in section 2. The data analysis proce- 
dure is described in section 3. Results, the catalogue of 373 
potential associations for 2FGL sources is given in section 
4 followed by concluding remarks. 



2 OBSERVATIONS 
2.1 Observing Program 

The areas of 95 per cent probability of 2FGL localisation, 
which corresponds to 2.448a", are typically 3-6 arcmin. They 
are still too large for a blind search with VLBI using tradi- 
tional approaches. Therefore, we organised observations in 
several stages. At the first stage we observed the sources 
with connected interferometers. Very Large Array (VLA) in 
the northern hemisphere, and ATCA in the southern hemi- 
sphere. Observations are made in remote wings of the wide- 
band receiving system that covers 4.5-10.0 GHz. The sky 
distribution of target sources is presented in Figure 4. 

We split the source list into two parts. A list of 411 
sources with declinations lower than +10° was observed with 
the ATCA. A list of 215 sources with declinations higher 
than 0° was observed with the VLA. This list includes 35 
overlapping sources in the declination band [0, +10°] and 18 
sources tentatively associated with supernova remnants and 
pulsar wind nebulae. Results of VLA observations will be 
reported in a separate paper (F. Schinzel et al., (2013) in 
preparation) . 

Observations at two frequencies allow us to determine 
spectral indices. Spectral index is defined as S* ~ where 
5* is the flux density, and / is the frequency. Sources with 
spectral index a greater than —0.5 will be considered as 
likely associations. Bright detected sources will be followed 
up with the VLB A and LB A. Detection of emission at par- 
sec scales will firmly associate them with the compact core 
regions of AGNs. 
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2.2 ATCA observations 

We observed a list of 411 target sources with ATCA on 
September 19-20 2012 for 29 iiours. We used automatic 
scheduling software sur_sked originally developed for VLBI 
survey experiments (Petrov et al. 2011b). In order to max- 
imise the observing efficiency, we reduced the correlator cy- 
cle time from the usual 10 second to 6 seconds, and ob- 
served in mosaic mode to minimise correlator overheads 
when changing sources. We tried to observe each target 
source in three scans 24 second long each separated by at 
least three hours in a sequence that minimises slewing time. 
In fact, we observed 80 sources in 4 scans, 239 sources in 
3 scans and 92 sources in 2 scans. At the beginning of 
the experiment we observed primary amplitude calibrator 
1934—638 for 20 minutes. The automatic scheduling process 
checked that for every target source a phase calibrator was 
inserted into the schedule that satisfy two criteria: the arc 
between the target and the calibrator should be less than 
20° and a calibrator should be observed within 20 minutes 
of a target. We used a pool of 1464 compact sources with 
5 < +20° with correlated flux densities > 200 mjy. The 
scheduling process picked up 207 sources. Many adjacent 
target scans reused the same phase calibrator. Each calibra- 
tor scan was observed for 18 seconds. 

The array was in the so-called II214 hybrid configu- 
ration with baselines in a range of 31-214 meters between 
the inner 5 antennas and ~4.4 km between CA06 and the 
inner antennas. We observed simultaneously in two bands, 
both 2GHz wide, centred on 5.5 and 9.0 GHz. In data re- 
duction the band edges were excluded and the ranges [4.58, 
6.42] GHz (thereafter 5 GHz band) and [8.09, 9.92] GHz 
(thereafter 9 GHz band) were used. The data were recorded 
at both polarizations with the Compact Array Broadband 
Backend (Wilson et al. 2011). 



3 DATA ANALYSIS 

The ATCA correlator provided us the output'^ with fre- 
quency resolution 1 MHz and time resolution 6 s. We used 
the software package Miriad (Sault, Teuben & Wright 1995) 
for development of the analysis pipeline. The data process- 
ing had the following steps: 

(a) We split the data into subbands and into scans. We 
discarded data from remote CA6 station for imaging analy- 
sis. 

(b) We analyzed observations of all the sources and 
flagged the data affected by radio interference using the task 
pgflag. 

(c) We used the Miriad task mfcal for computing band- 
pass calibration using the data of 1934—638 with solution 
interval 2 minutes. This bandpass was applied to all the 
data. 

(d) We determined antenna gains and polarization leak- 
age of phase calibrators using the Miriad task gpcal and we 
also corrected fiux scale with the task gpboot. Then these 
complex gain factors were applied to target sources. 



Available at http://astrogeo.org/vl/aofusrpfits 
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Figure 5. Example of a 9 GHz image of around 2FGL source 
J1353.5— 6640. Flux densities were not corrected for beam atten- 
uation. Objects 3-8 were flagged out. The gray area in the center 
of the image shows the error ellipse with the probability 95 per 
cent to find there J1353.5— 6640. 



(e) All calibrated visibilities of a target source were 
merged into one file. The visibilities were inverted using the 
Miriad task invert with the multi-frequency synthesis algo- 
rithm. Typical image size was 1024 x 1024 with a pixel size of 
2.4 arcsec although for some target sources we changed the 
pixel size later during data analysis. We tapered the data ac- 
cording to Brigg's visibility weighting robustness parameter 
factor 0.5 during inversion. 

(f) We cleaned the data using the Miriad task mfclean. 
The default CLEAN gain was set to 0.01 and the default 
maximum number of iterations was set to 192. The default 
clean region was 93 per cent of the image. 

(g) We searched for point sources in a clean image using 
the task imsad with the flux densities greater than 5.5 image 
noise root mean square (rms). 

(h) We performed the non-liner least square fit for po- 
sitions and flux densities of detected sources to calibrated 
visibilities using the Miriad task uvfit. We restricted the fit- 
ting parameter space to the vicinity of parameters found 
by imsad. This procedure also provided us uncertainties of 
parameters. We ran this procedure for each subband sepa- 
rately. 

At this point we checked visually every clean image. 
We have developed a web-application that allows us to in- 
teractively examine each image and flag those point sources 
found with program imsad that either seem spurious or lie 
far away from the pointing direction. For 7 per cent images 
we manually adjusted parameters for inversion and clean- 
ing, such as image size, pixel size, offset of image center, 
and the number of CLEAN iteration. Figure 5 provides an 
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illustration'' . Sources 4-8 have a signal-to-noise ratio (SNR) 
in the range 5.5-6.7 — not sufficient for a reliable detec- 
tion in the absence of other information. There are two 
sources within the 95 per cent localisation error ellipse of 
J1353.5— 6640. Source 1 has spectral index —0.46, while 
source 2 has spectral index —1.37. There is a strong X-ray 
source IRXS J135341. 1—664002 within 7 arcsec of source 1. 
Tsarevsky et al. (2005) observed with VLT its optical coun- 
terpart within 1.7 arcsec of source 1 and reported a feature- 
less spectrum. He suggests its classification as a BL Lac type. 
This additional information presents a strong argumentation 
in favor of association of source 1 with 2FGL J1353.5— 6640. 

Source 3 was flagged out because it lies 10.3 arcmin from 
the pointing direction. The beam power at 9 GHz drops by 
a factor of 2 at 2.7 arcmin and has the first minimum at 
6.5 arcmin. In fact, source 3 lies within 1 arcsec, i.e. within 
its position error, from J1353-6630 detected in AT20G sur- 
vey with ATCA at 20 GHz (Murphy et al. 2010). That cat- 
alogue reports its flux density at 8 GHz 322 mjy. Flux den- 
sity from our observations appeared 77 times less because 
it was detected by a sidelobe of the beam. We did not try 
to calibrate source flux densities for objects that had offsets 
from the pointing direction beyond where the beam power 
falls below 0.2: 400 arcsec at 5 GHz 250 arcsec at 9 GHz. 
For sources that lay within these limits we corrected their 
flux densities for beam attenuation P{9) using approxima- 
tion of the beam pattern measured by Wieringa & Kesteven 
(1992) in the form of 

(4) 

k = 

where 9 is the distance from the pointing direction, / is the 
observing frequency, and Ck are coefficients that approxi- 
mate beam pattern measurements. 



4 RESULTS 

In summary, we found 626 objects in 267 fields. Of them, 
210 objects have position offsets with respect to 7-ray 2FGL 
objects exceeding 2.448(t, i.e. beyond the areas of 95 per cent 
probability of their localisation, and they were discarded. 

The rms of images is typically in a range of 0.15- 
0.25 mJy. The FWHM size of the restored beam is typi- 
cally 35 arcesec at 5 GHz and 20 arcsec at 9 GHz. The 
detection limit is 1.8 mJy for sources in the center of the 
field of view and 9 mJy at the edge of the field of view, 
at 6.5 arcmin. There are 23 fields, for instance the field 
around 2FGL source J1619.7— 5040, that have a number of 
resolved objects. The field around J1619.7— 5040 lies in the 
HII region G332. 8-0.6 (Kuchar & Clark 1997). Although 
the data analysis pipeline selected sources that it consid- 
ers as "points", these are rather hot spots in an extended 
Galactic object than separate compact extragalactic sources. 
Imaging of an extended object using 2-4 scans with a five- 
element array gives inconclusive results. We flagged sources 

All 822 images derived from analysis of these observations are 
available online at http://astrogeo.org/aofus 



from such fields as possibly extended. Extra care must be 
taken when dealing with these objects. 

We cross-matched all remaining sources against NVSS 
(Condon et al. 1998), SUMSS version 2.1 of 2012 February 
16 (Bock, Large & Sadler 2003; Mauch et al. 2003), MGPS- 
2 (Murphy et al. 2007), and WISE (Wright et al. 2010) cata- 
logues. NVSS catalogue is derived from VLA observations at 

1.4 GHz, SUMSS and MGPS-2 are derived from processing 
observations at 0.843 GHz with the Molonglo Observatory 
Synthesis Telescope. These catalogues have similar resolu- 
tions ~40 arcsec and are complimentary to each other. We 
used 20 arcsec search radius for matching with an NVSS ob- 
ject and 2(T position uncertainties of its ATCA coordinates 
for matching a WISE object. 

We have 136 matches against WISE, 174 matches 
against NVSS, 57 matches against SUMSS, and 48 matches 
against MGPS-2. Since the average position error of WISE 
is rather small, 0.25 arcsec, the mathematical expectation of 
the number of background WISE sources that fall within 2a 
error ellipse of ATCA position in a case if the ATCA detec- 
tions and WISE objects are physically unrelated is only 10. 
We found matches for 13 times more sources, roughly 1/3 of 
our list. Thus, we conclude that the majority of ATCA and 
WISE matches are the same sources. 

The NVSS catalogue has objects with declinations > 
-40°. We have 227 detections with 5 > -40°, i.e. 77 per 
cent of our objects have a counterpart in NVSS. Among 190 
detected objects with declinations < —40°, 55 per cent have 
counterparts in SUMSS and MGPS-2. We have less coun- 
terparts with declinations below —40° because these cata- 
logues are not as deep as NVSS: their limiting peak bright- 
ness is 6 mJy beam^^ at 5 ^ —50° and 10 mJy beam~^ 
at (5 > —50°, while NVSS has the limiting peak brightness 

2.5 mJy. 

We used these matches to re-calibrate position errors. 
We selected 92 matches with flux densities at 1.4 GHz 
brighter than 15 mJy. NVSS positions of such sources are 
known with errors less than 1 arcsec. We formed the dataset 
of position differences normalised to their standard devia- 
tions which is a sum in quadrature of NVSS position un- 
certainties and position uncertainties from our ATCA ob- 
servations. We sought the reweighting model in a form 
y^(r"ap"+"y^ to bring the the averaged normalised resid- 
uals close to one. We found the following scale-factors and 
error floors: 1.45 and 0.5 arcsec for right ascension and 1.80 
and 0.75 arcsec for declination. We inflated position errors 
according to this reweighting model. 

For sources detected at both subbands independently, 
we determined spectral indices from their flux densities. A 
number of sources was detected at the low subband only. 
For strong sources we were able to determine their spectral 
index within the band [4.58, 6.41] GHz by splitting the sub- 
band into 8 segments and re-running the last step of the 
analysis procedure for each segment: least square flt for po- 
sitions and flux densities in the vicinity of detected sources. 
Then we fitted the straight line into the dependence of fiux 
density S on frequency / as S{f) = -I- k(/ — /o) where 
fo is the weighted mean frequency. The spectral index a is 
determined as q « k/o/Sq. 

In general, spectral index from both subbands and from 
the low subband spectra is not necessarily coincide, because 
systematic errors affect the estimates in a different way and 
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Table 1. The first 8 rows of 142 objects detected at both 5 and 9.0 GHz sub-bands. The table columns are explained in the text. The 
full table is available in the electronic attachment. 
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Figure 6. Comparison of source spectral indices from dual sub- 
band observations (horizontal axis) with the spectra determine 
from the low sub- band only (vertical axis). 



some sources may have spectra that deviate from the power 
law. In order to compare a consistency of spectral index es- 
timates, we computed spectral indices using the low sub- 
band spectra for the sources that were detected at both 
bands as well. The scatter plot of this comparison is shown 
in Figure 6. Results of this comparison demonstrate a rea- 
sonable consistency of spectral indices estimates using both 
approaches. 

In order to evaluate the significance of associations, we 
computed the likelihood ratios A defined as a ratio of the 
probability that a radio counterpart will be found inside the 
disk of radius d equal of their separation to the probability 
to find a background radio with flux S or greater outside the 
same disk: 



A = 



iV(l)S"'d74 



(5) 



where n is normalied arc d between radio and 2FGL po- 
sitions, A^(l),p are parameters of log A'^ (log S) dependence 
that describes the number of sources A'^ with radio flux 
densities greater than 5. Among large radio surveys, GB6 
(Gregory et al. 2003) catalogue made at 4.85 GHz with 
FWHM of the beam 3 arcmin, is the closest match to our 
ATCA observations at 5 GHz. We built a cumulative his- 
togram log N-log S over 75,162 GB6 sources at 6.06 sr and 
fit a straight line into it in the interval of flux densities [30, 
810] mJy. We got N(Uy) = 965 and p = -1.3475 from this 
analysis. We should note that these parameters are different 
than those we derived from analysis of VLBI results, because 
they are related to a different population of radio sources. 



The likelihood ratio describes our knowledge of how likely 
a radio source with a given flux density can be found at a 
given distance by chance. Sources with large A are less likely 
located close to the 7-ray counterpart by chance. 

The list of 416 sources falls into three categories: 

(1) Category I: 142 objects detected at both 5 GHz and 
9.0 GHz subbands within 2.7 arcmin of the pointing direc- 
tion. We provide in Table 1 7-ray source name; lAU name 
of detected radio source; tentative association status, esti- 
mates of its J2000 coordinates followed by la uncertainties 
in arcseconds (uncertainty of right ascension is not scaled by 
cos S factor), Aa and AiS; flux densities at 5 GHz and 9 GHz 
in mJy, F5 and Fg, corrected for beam attenuation followed 
by their standard deviations, A5 and Ag; and spectral index 
and its error, Sp and ASp. We provide distance of a source 
from the pointing direction D, followed by Na, the ratio of 
this distance to its standard deviation derived from reported 
2FGL position localisation errors. If the source was associ- 
ated with an object either from NVSS, SUMSS or MGPS-2 
catalogues, its flux density at 1.4 GHz (NVSS) or 0.843 GHz 
(SUMSS and MGPS-2) is shown in column Fi followed by 
the 1 GHz catalogue code {N for NVSS, S for SUMSS, M 
for MGPS-2), followed by the source identifler in that cat- 
alogue. If the source was associated with WISE object, its 
WISE source ID is shown in the last column. 

Golumn 3 shows two flags: "a" if a source has likelihood 
ratio greater than 10, and therefore, considered a likely as- 
sociation, "e" if a source is extended, and "f" if a source has 
the spectral index flatter than —0.5. 

(2) Category II: 210 objects detected only at 5 GHz 
subband within 6.5 arcmin of the pointing direction. We pro- 
vide in Table 2 estimates of fiux density at 5 GHz corrected 
for beam attenuation. The contents of Table 2 is similar to 
Table 1, except columns Fg and Ag and which are missing. 
The spectral index is computed over the sub-band [4.58, 
6.41] GHz. The spectral index estimates with uncertainties 
greater than 0.3 are omitted. 

(3) Category III: 64 objects either detected beyond 
6.5 arcmin of the pointing direction or detected only at 
9 GHz band. Since calibration for beam attenuation becomes 
uncertain at large distances, we can only provide a low esti- 
mate of their flux density: 10 mJy. Table 3 lists these sources. 
Its contents similar to Table 2, except columns F5 and A5 
which are missing. 

As can be seen from Figure 7, many of the unassoci- 
ated Fermi sources are at low Galactic latitude. For the 
30% of the sources in Tables 1-3 which lie within 1 degree 
of the Galactic plane, contamination of our ATCA images 
by Galactic HII regions, supernova remnants and planetary 
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Table 2. The first 8 rows of 210 objects detected at 5 GHz subband only. The table columns are explained in the text. The full table is 
available in the electronic attachment. 
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Table 3. The first 8 rows of 64 objects detected at sidelobes of the 5 GHz beam, or detected only at 9 GHz. The table columns are 
explained in the text. The full table is available in the electronic attachment. 



2FGL name 


lAU name 


Fl 


Q 


5 


Aa 


A5 


D 


No- 


Fx 


Ci 


IGHz ID 


WISE ID 










hr min sec 


O / // 


II 


II 






mJy 








(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


(12) 


(13) 


J0414.9- 


0855 


J0415-0854 




04 15 23.27 


-08 54 26.5 


3.6 


3.5 


6.9 


2.1 


10.3 


N 


041523-085425 


J041523. 13-085424.4 


J0540.1- 


7554 


J0539-7601 




05 39 18.29 


-76 01 30.6 


13.7 


4.4 


8.0 


2.2 


79.0 


S 


053916-760131 


J053916. 48-760130.3 


J0540.1- 


7554 


J0540-7601 




05 40 10.51 


-76 01 54.0 


11.8 


3.9 


7.7 


2.2 


37.0 


S 


054011-760156 


J054011. 59-760154.3 


J0540.1- 


7554 


J0541-7601 




05 41 22.48 


-76 01 09.6 


18.9 


5.8 


8.1 


2.4 


31.0 


s 


054123-760108 


J054123. 16-760107.6 


J0547.5- 


0141 


J0547-0133 




05 47 20.86 


-01 33 30.9 


1.4 


1.8 


8.1 


1.4 


22.1 


N 


054720-013327 


J054720.85-013329.9 


J0555.9- 


4348 


J0555-4346 




05 55 16.04 


-43 46 30.7 


1.8 


1.4 


7.1 


2.2 


144.0 


s 


055516-434631 




J0555.9- 


4348 


J0555-4345 




05 55 19.57 


-43 45 28.9 


2.2 


1.7 


6.8 


2.1 


101.0 


s 


055519-434534 




J0600.8- 


1949 


J0600-1950 


a 


06 00 59.64 


-19 50 39.4 


1.0 


1.0 


2.0 


0.5 


96.7 


N 


060100-195049 





OETT— 1 r- 

o 

o : 




-60 -40 -20 20 40 60 

Galactic latitude (deg) 

Figure 7. Dependence of flux density of detected sources of cat- 
egories I and II on galactic latitude. Flat spectrum sources are 
shown with green hollow circles. Other sources are shown with 
filled blue discs. 



nebulae is common so careful -work -will be needed in the final 
identification process. As expected, flat-spectrum sources do 
not exhibit concentration to tlie Galactic plane. 
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5 CONCLUDING REMARKS 

We found 416 counterparts within 95 per cent probability of 
localisation of 267 2FGL sources (Tables 1-3). For 69 2FGL 
sources more than one counterpart was found. Among 416 
counterparts, 95 have a spectrum flatter than —0.5, 75 have 
a spectrum steeper than —0.5, and for 246 objects no spec- 
tral index has been determined. The flat spectrum sources 
are considered as tentatively associated with Fermi objects. 
For 107 sources the probability of association is more than 
10 times exceeds the probability to find a background source 
with that flux density by chance. These flags help to iden- 
tify a subset of more probable associations. However, taken 
alone, they do not establish a firm association. For instance, 
object J1045— 5941 has spectral index 1.07, but it is a Galac- 
tic object (77 Carinae), not a blazar. Some compact HII re- 
gions and planetary nebulae may have flat spectrum. Ad- 
ditional information, f.e., from observations at other wave- 
lengths will help to identify a source associated with a 7-ray 
object. 

Similarly, steep spectrum does not rule out extragalac- 
tic nature of an object. Kovalev (2009) presented several 
examples of steep-spectrum sources associated with Fermi 
objects. Detection of emission from parsec-scale will allow 
us to determine the nature of a radio source. 

There are 46 objects from the fields that that have a 
number of resolved objects. Many of these fields are asso- 
ciated with HII regions. Objects from these fields may be 
hot spots of extended objects, not separate sources. Only 4 
of them, i.e. 9 per cent, are associated with WISE objects, 
while 35 per cent remaining sources have counterparts from 
WISE catalogue. 

We did not find radio counterparts for 144 or our tar- 
gets, roughly 1/3 of our list. Since for many sources the 
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2.448(7 2FGL position error ellipse exceeded the field of view 
of 6.5 arcmin at 5 GHz, we cannot claim there is no source 
brighter than 10 mjy in the field. 

In the second phase of the project we are going to re- 
observe these fields in the mosaic mode. This will allow us to 
detect all the sources in the 95 per cent probability of 2FGL 
localisation that are brighter than 2 mJy. We are also going 
to re-observe with ATCA those sources that were detected 
at 5 GHz, but are not detected at 9 GHz, and had posi- 
tion offsets exceeding FWHM of the 9 GHz beam. We will 
re-observe sources listed in Table 3 with pointing to their 
positions found from this survey as well. 

In the third phase of the project we will observe de- 
tected sources with LBA and determine correlated flux den- 
sities from regions smaller than 50 mJy. This will allow us 
to associate detected radio sources with blazars. 
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